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I. INTRODUCTION 



Pions are abundantly produced in heavy-ion collisions at beam energies of 1 A GeV and 
above. The main pion production channel is via the resonance excitation in nucleon-nucleon 
collisions: NN — > NR, R — > Nit, where R stands for A(1232) or higher baryonic resonances. 
As well established in nuclear transport theory [[[], EL pions propagating through nuclear 
matter experience the chain of absorption-production processes tciNi — > — > 7r 2 A^ 2 ~~ * 
R2 etc due to the large resonance irN scattering cross section. In other words, a 

pion with momentum ~ 0.3 GeV/c acquires a width in nuclear matter. The pion spectral 
function in nuclear matter is, therefore, different from the one in vacuum. Actual calculations 
within the (generalized) A-hole model || |j, ||, |, indeed have shown that the pion spectral 
function is quite broad. 

The off-shell nucleon effects in heavy particle (N*(1535)) production have been studied in 
||, where it has been shown, that the ground state correlations in colliding nuclei enhance the 
near-threshold production substantially. Beyond the ground state correlations, the influence 
of the off-shell nucleon propagation on heavy particle production has been found to be 
unimportant in || , while recent BUU calculations |], [TO] have demonstrated an appreciable 
enhancement of the light meson (n,K + ) production by off-shell nucleons. 

The transport simulations with in-medium pion and A properties have been recently 
performed by Helgesson and Randrup in a series of works (see and Refs. therein) 
employing the quasiparticle approximation for pions. The soft pion enhancement due to the 
in-medium pion modifications has been observed in |TT]] in agreement with BUU studies in 



12] , [13 ] . As has been pointed out in |TTJ, the overall effect of the in-medium modifications 



on the pion observables is small due to the fact, that the observed pions are emitted from 
the surface region of the system, where the in- medium modifications are not important. 
However, by multistep processes, pions influence other particles, e.g. kaons. Thus, the in- 
medium pion modifications could manifest themselves in the changes of the yields of the 
other particles. 

The present work is the first attempt to include the pion spectral function into the BUU 
transport theory. The structure of the paper is as follows. In Sect. II we derive the pion 
spectral function in nuclear matter at finite temperature on the basis of the A-hole model. 
Sect. Ill contains the equation of motion for off-shell pion propagation. In Sect. IV we 
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present the results of the BUU calculations for the pion production in Au+Au collisions at 
1 A GeV. The summary of the work is given in Sect. V. 



II. THE PION SPECTRAL FUNCTION IN NUCLEAR MATTER 

In the A-hole model the pion polarization function H(k) is given by the loop shown in 
Fig. |l], where we neglected the exchange it-diagram. For the evaluation of IT(/c) we use a 
simple Lagrangian of the form: 

ZttNA = ^^V^T^tv^TT + h.C. , (1) 

where T is the isospin transition operator (1/2^3/2) [|J], if>^ is the Rarita-Schwinger field 
of the A-resonance, ipN and 7r are the nucleon and pion fields, m n = 0.138 GeV is the pion 
mass. The coupling constant f n NA (— 2) will not enter explicitly in final expressions for 
H(k) (see Eqs. (jI7|),(|I1f),(f2"0])), therefore, we do not need to specify its value. Using Feynman 
rules and the Lagrangian the pion polarization function in the case of isospin symmetric 
nuclear matter may be written as: 

f d*p 



- ill(k) = / (f^ASp (GT(P + k)G (p)) , (2) 



G (p) = (p' + m N ) (- \— - + 2m5(p 2 - m 2 N )9(p )n(p)) (3) 



where 



is the nucleon propagator in noninteracting nuclear matter (c.f. 

n(p) = exp((E N (p)-n)/T) + l (4) 
is the occupation number written for simlicity in the rest frame (r.f.) of nuclear matter. 



E N (p) = ^pi + m% (5) 

is the free nucleon energy with = 0.938 GeV being the nucleon mass. For later use we 
define here also the free A and pion energies: 



E A (p) = Jp 2 + ml , (6) 



E n (k) = JW + ml (7) 
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with the pole delta mass ira = 1.232 GeV. The Rarita-Schwinger A-resonance propagator 
(*a(p) * s n °t unambiguously defined in the literature. We will use an expression from Ref. 

aW p 2 - mi + iv^r A (p2) ^ 3 7 7 3p2 + 3v ^ J' W 

which is similar to the one from Ref. |J4| with an addition of the term proportional to the 
total width of the A-resonance T&(p 2 ) in the denominator and with a replacement of the 
pole delta mass m A by y/p 1 in the numerator and in the tensor part. 

The nucleon propagator (Q) contains the usual vacuum part oc (p 2 — m 2 N + ie)~ l and the 
in-medium part oc 5(p 2 — m 2 N ). The latter reflects rescattering of a pion on real nucleons. 
The vacuum part contributes to ImlT(fc) only for k 2 > (m^ + m A ) 2 , which corresponds 
to the decay of an off-shell pion to the A-resonance and antinucleon. Obviously, the pions 
produced in the baryonic resonance decays can not match this condition. Therefore, we will 
keep only the in- medium part of the nucleon propagator. Substituting Eqs. into Eq. 

(0) and keeping only the density dependent (oc n(p)) term in the nucleon propagator we get 
the following expression for the pion polarization function: 

where s = (p + k) 2 and the nucleon is on its mass shell, i.e. p 2 = m 2 N . In the nonrelativistic 
limit |k| CniAf, |p| <C and for small A off-shellness \y/s—m&\ -C m A Eq. (^) transforms 
to well known formula (c.f. Ref. 0): 

m ~ |k| 7 Wf n{v) E N { v ) + kP - e!( p + k) + ir A ( s ) • (10) 

In the spirit of the transport theory, however, one would like to express the pion polar- 
ization function in terms of the ir + p — > A ++ cross section 

< u*\ - 4tt 2sT a ^ wN (s, k 2 )T A (s) 
a ^ A++[S ' k } ~ q 2 (s,k 2 ) (s - m 2 A ) 2 + sT 2 A (s) ' W 

generalized here for the case of the off-shell pion with four-momentum k. 



q(s, k 2 ) = ^/(kp) 2 - k 2 m 2 N /V~s (12) 

is the cm. momentum of the pion and nucleon while Y^ 7r ^{s, k 2 ) is the decay width of 
the A-resonance. The decay width can, in principle, be easily calculated within the same 



Lagrangian of Eq.([l|) which produces the expression (c.f. Ref. fl4|): 

r£L^(*,* a ) = fl ^[ s : k P (W 2 - ™%k*M P + k) + m N VS) 



Vl'KVfv^s 2 



<Jq2(s,k 2 )+m 2 N + m N ) . (13) 



\2nm1 y/s 

However, we will use Eq.(|i~3|) only for the derivation of the pion polarization function in terms 
of the cross section. In numerical calculations of the pion spectral function we substitute in 
Eq.flTTD the off-shell extension of the vacuum width parameterization from Ref. [16] 



rw(*, k)-r ^-j —^ + q2 ^ k2) (14) 

with g = Qi'm'A, m n), To = 0.118 GeV and (3q = 0.2 GeV. The total width of the A-resonance 
in nuclear medium is taken as 

T A (s) = T A _^ N (s,ml) . (15) 
Substituting Eqs.([TT|),([T3|) into Eq.(|9|) one gets the expression: 

U{k) = 3 J (2nrE N ( P ) n{p) fZ£) ^A ++ («, fc ) • (16) 

Separating now real and imaginary parts we have: 

f 2d?p m N ,prei(s - m\) 4 2 
Ren = / ra(p) »^+ P ^A++(s,fc , (17) 

7 (27r) J £W(p) r A (s)Vs 3 

T tt/m f 2d 3 p m N 4 2 

Imll(fc) = - / /«-) \3 T7i / M p)p re i-a 7r+p ^ A ++(s,k ) , (18) 
J (Ait) 6 h N {p) 3 

where p re \ = qy/s/rriN is the pion momentum in the r.f. of a nucleon. By noting that 
4/3<7 7r + p _ > A++ = cTttp^a + Cnn->A we see that Eq.fllTf) can be simply rewritten in a more 
compact and intuitive form: 

ImTl(k) = -~( m f r p re ;(cr 7r p_ tA + 0"™-^)) , (19) 

where the averaging is performed with respect to the nucleon momentum and p is the nuclear 
matter density. Notice, that Eqs.(|T7|),(|T8|) are explicitly relativistically covariant, since the 
self-energy of a scalar particle is a Lorentz-scalar. 

As pointed out in the beginning of this section, in the derivation of the pion polarization 
function we have missed the w-diagram. The particular form (|8]) of the spin-3/2 propagator 



does not permit to reliably calculate the w-diagram contribution, since the tensor part has a 
pole at p 2 = (c.f. [113] for details). However, using the usual Rarita-Schwinger propagator 



(c.f. [|14|]) we have checked that the w-diagram contributes mostly in the space-like region 
k 2 < 0, which is not important for our study, since in BUU only the time-like pions are 
propagated. 

Finally, we correct the pion polarization function for the repulsive interaction of holes 
and A- resonances at short distances (c.f. ||): 

n ^'» - w^m (20) 

with g' = 0.5 being the Migdal parameter. Here the pion four-momentum k = (k°, k) has 
to be taken in the r.f. of nuclear matter. In that frame, due to the isotropy of the nucleon 
Fermi distribution, the polarization function n c depends only on the pion energy and on the 
absolute value of the pion momentum H c (k) = U c (k°, |k|). 

Fig. H shows the real and imaginary parts of the pion optical potential 

UW) - , (21) 

where Sk is the solution of the dispersion relation (c.f. J7J) 

4 = i^(k)+Ren c (£ fe ,|k|) . (22) 

We see that the real part of the optical potential changes sign and the absolute value of the 
imaginary part reaches its maximum practically at the same momentum |k| ~ 0.380 GeV/c 
for normal nuclear matter case: p = po, T = 5 MeV (solid lines), where po = 0.16 fm -3 is the 
nuclear saturation density. This can be qualitatively understood from Eqs.flT7|),(pif) since, 
neglecting for simplicity the Fermi motion, Ren = and ova^a^) is maximal at s = m\. 
The estimated pion momentum is |k| = 0.300 GeV/c, which is somewhat different from 
the numerical value above due to the Fermi motion and the Migdal correction. Increasing 
temperature (dashed lines) smeares out the jump of Re U and decreases the peak value of 
— Im U while shifting the peak to lower values of the pion momentum. Increasing density 
(dot-dashed lines) acts in the opposite direction. These findigs are in an overall agreement 
with the calculations of Ref. (see Fig. 2 of this Ref.). An essential difference is visible 
only at high momenta |k| > 0.4 GeV/c where the potentials of Ref. are more repulsive. 
The pion optical potential is, however, sensitive to the choice of the nucleon and delta 



propagators (relativistic in our work and nonrelativistic in Ref. [0]), to the presence of the 
^-diagram contribution missed in our calculations for simplicity and - most important - to 
the choice of the A-resonance width. In this exploratory work we do not intend to study in 
detail all these effects. Fig. ^| has a purpose to demonstrate only that our simple spectral 
function calculation is in the reasonable agreement with a state-of-art finite temperature 
calculation of Ref. [|7|. 

Next we determine the pion spectral function 

1 1 ImIl c (A;) 

An{k) = --l m G(k) = -- (p - ^ _ Renc( ^ ))2 - (Im n c (A;))2 ( 23 ) 

with G(k) = [k 2 — m 2 — I1 C (A;)) _1 being the pion propagator. The exact spectral function 
must satisfy the normalization condition 

d(k°) 2 A n (k) = 1 . (24) 



o 

We checked that the spectral function ([23|) satisfies the condition ([24]) with the accuracy of 
10 %. 

Fig. ^ shows the spectral function vs pion invariant mass squared M 2 = k 2 and momentum 
|k| for normal nuclear matter (density p = po and temperature T = 5 MeV). There are two 
peaks of A as a function of M 2 at fixed momentum |k|. At |k| < 0.4 GeV/c the lower M 2 
sharp peak corresponds to the pion-like collective mode, while the higher M 2 broad peak 
reflects the A-hole mode. With increasing momentum the lower M 2 peak becomes wider 
and enters the negative M 2 (space-like) region, while the higher M 2 peak gets narrower and 
asymptotically approaches the free pion limit (M 2 = m 2 ~ 0.02 GeV 2 ). Therefore, at high 
momenta |k| > 0.4 GeV/c the higher M 2 peak corresponds to the pion-like mode and the 
lower M 2 one describes the A-hole excitation. This level crossing is better illustrated in 
Fig. |] where along with the peak positions of the spectral function in the (|k|, a;)-plane (full 
circles) we present also the free pion dispersion relation uo = E n (li) (solid line), the A-hole 
mode 



E N A(k) = yJu A {k){u A {k) + g'C) (25) 

with co>A(k) = E A (k) — rriN and C = 8/9(/ 7r ATA/^ 7 r) 2 p (c.f @) (dashed line), and the solution 
Ek of the dispersion relation ( |22"D (dash-dotted line). The peak positions close to the free pion 
solution are very well described by the dispersion relation (E3), which has a characteristic 
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jump at |k| = 0.38 GeV/c. We also observe, that broad peaks at higher (lower) energy for 
small (large) momenta correspond indeed to the A-hole excitation. 

The high-momentum structure of the pion spectral function (|k| > 0.4 GeV/c) would 
be modified by taking into account the higher resonances in the pion polarization function. 
Therefore, the pion peak at high momenta is expected to be smeared out by the higher 
resonances. Nevertheless, as shown in Ref. JT7| for the central heavy-ion collisions at 
energies about 1 A GeV, the higher resonance contribution dominates in the pion yield 
at the cm. pion kinetic energy above 0.6 GeV which corresponds to pion cm. momentum 
above 0.7 GeV/c. Thus, we expect that the inclusion of higher baryonic resonances in the 
pion spectral function will not strongly modify the results of our work (c.f. Fig. [7] later). 

In order to estimate an effect of the in-medium modifications of the A-resonance width, 
we replaced Ta(s) in Eqs. by Ta(s)(1 — n) +r sp , where T sp = — 21m W sp = SOp/po 

[MeV] is the phenomenological spreading width given by the imaginary part of the A- 



spreading potential ||18|| . n is the angle averaged nucleon occupation number. The A- 
spreading potential takes into account the quasielastic scattering, two- and three-body ab- 
sorption contributions |l9, 21J. In Fig. || we compare the spectral functions at p = 2p and 
T = 70 MeV calculated with vacuum (dashed lines) and in-medium (solid lines) A-resonance 
width. We see that the A-spreading potential smeares-out the peaks of the spectral function. 
There is a tendency of the strength concentration near M 2 = m 2 = 0.02 GeV 2 . Correspond- 
ingly, the strength at negative and large positive M 2 is reduced somewhat by in-medium 
A-width. 



III. KINETIC EQUATION FOR OFF-SHELL PIONS 

For the description of the off-shell particle dynamics we will use the spectral distribution 
functions (c.f. 0) F Q (r, t, p, M 2 ) with a being the particle type (a = N, R, 7r for nucleons, 
resonances and pions respectively). The number of particles of type a in the phase space 
element d 3 rd 3 pdM 2 / (2tt) 3 is F a d 3 rd 3 pdM 2 / (2ir) 3 . The following relation holds: 

F a (r, t, p, M 2 ) = / a (r, t, p, M 2 )A Q (r, t, p, M 2 ) (26) 

where / a (r, t, p, M 2 ) is the usual phase space density, and A a (r, t, p, M 2 ) is the spectral 
function. The pion spectral function has been defined earlier by Eq . (|23|) . The spectral 



S 



function of a baryonic resonance is taken as the vacuum one: 



A (M 2 ) - 1 MT R {M 2 ) 
R[ } ~ n(M 2 -m 2 R ) 2 + M 2 T 2 R (M 2 ) [Z<) 



with m R being the pole mass and T R (M 2 ) being the vacuum width of the resonance [|16 



The effect of the nucleon off-shellness on the pion and kaon production has been studied in 
Ref. |J. In the present work we treat, however, nucleons on-shell in order to underline the 
genuine pion off-shellness effects on the observables. Thus, we use 

A N (M 2 ) = 5(M 2 - m%) (28) 

and everywhere below the nucleon four-momentum is pn = (En(pn), Pn)- 

It is instructive to write down explicitly the coupled transport equations for resonances 
and pions: 

( d dH* f d dH* d \ , 

(a + -^S--* £ flSj^ r ' f ' P *^ - (29) 

^ + ^jF 7r (r,t,k,M|) + L 0// _ s ^[F 7r ] = I n , (30) 

where H^(r, t, p R ) is the mean field Hamiltonian acting on resonances: in this work we take 
it the same as for nucleons. I R and 1^ are the collision integrals. The term L ff_ s h e u[F n ] 
on the l.h.s. of (|30"D, which will be specified later on, is introduced in order to bring pions 
back on the mass shell once they reach the vacuum. The collision integrals are expressed in 
terms of gain (<) and loss (>) terms as follows: 



,< _ ,> = v- J_ f 9Nd 3 p N r d 3 kdM 2 
R R ~f^2p%J (2n)*2E N (p N ) J (2tt)32A;0 



x (2ir) 4 5^(p R -k-p N )\M R - 



•Ntt\ 



|2 



x {F w (k, M 2 )f N (p N )A R (M 2 R ) (1 - f R (p R , Ml)) 

- F R (p R , M^A^h, M 2 )(l + U(k, M 2 n ))(l - f N (p N ))} , (31) 



r < r > _ 1 /" 9Rd 3 PndM R r g N d 3 p N 



* w ** ^ 5 J (2n)*2jP R J (2itf2E N (3> N ) 

x (2ix) A 5^\p R -k-p N )\M R ^\ 2 

x {F R (p R ,M 2 R )A n (k,M 2 )(l + U(k,M 2 ))(l-f N ( PN )) 

- F n (k,M 2 )f N (p N )A R (M 2 R )(l - f R {p R ,M 2 R ))} , (32) 

where we dropped space and time variables for brevity. In Eqs. ([[]]), I-^-b^NttI 2 is the 
invariant matrix element squared averaged over spins of initial and final particles, g a = 



2s a + 1 (a = N,R) is the spin degeneracy factor. For simplicity, we did not display the 
contributions of NN <-» NR processes to I R and of the direct processes NN «-> NNir to 7^. 
However, these processes are always included in BUU. Sums in Eqs.(|31|),(|3"2"|) run over all 
possible isospin states of the particles. Noting that the resonance decay width to a nucleon 
and an off-shell pion of mass M n in the resonance r.f. is given by the expression 

1 f gNd 3 pN f d 3 k 



R 



2M R J (2n) 3 2E N (p N ) J (2n) 3 2k° 

x (2n) 4 6^(p R - k - p N )\M R ^ N7T \ 2 = I <m c . m 1 \M R ^ Nw \ 2 qg N (33) 

J 62tt z M r 

and using the fact that the spin-averaged matrix element does not depend on the center-of- 
mass scattering angles, one can rewrite Eq.([H]) as follows: 

W ' dM 2 ^T R ^(M 2 R ,M 2 ) 
" T J Pr 

^{^(k, M 2 )f N (p N )A R (M 2 R )(l - f R (p R , Ml)) 

- F Jl (p Jlj Jl^A r (k J ^)(l + / )r (k ) A^))(l-/ JV (p JV ))} . (34) 

Introducing the cross section 



IR - 

TV,- 



A 2 

a^ R (M 2 R , Ml) = -^IW^MS, M 2 )M R A R (M 2 R )^ (35) 
Q 9n 



the pion collision integral can be also simplified as 

g N d 3 p 

R.N 



\- f 9Nd 3 pN ,,, 2 
2^ J ( 2 7t) 3 v *nVkN->r{M r , MJ 



(27T) 

x {f R (p R , M 2 R )A n (k, M 2 )(l + U(k, M 2 ))(l - f N (p N )) 

- F^,Ml)f N (j> N ){l-f R {p Rl M 2 R ))} (36) 



with ivtv = y (kp N ) 2 — M 2 m 2 N /(k°E N (p N )) being the relative velocity of the pion and 
nucleon. 

In the following we neglect for simplicity the second-order effect of the pion off-shellness 
on the resonance width. Thus, we replace in Eq.(^) and, to maintain detailed balance, also 
in Eq.(j35|) 

T R ^ Nn (M 2 R , Ml) - T R ^ Nn (M 2 R , ml) , (37) 

i.e. the off-shellness of the pion in the resonance decay vertex is neglected. We have to 
comment that this simplification has been done only in the BUU implementation. (In the 
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pion polarization function calculation the full off-shell decay width (|i~4"D has been applied.) 
Correspondingly, the mass squared of the outgoing pion in the resonance decay has been 
sampled in the interval [0; — m^) 2 ] according to the distribution oc ^(k, Ml). In the 
case of the direct (background) pion production NN — > NNtt the value of Ml has been 
sampled in the interval [0; (\^s — 2m^) 2 } using the same distribution. 

The propagation of the off-shell pion between its production in the resonance decay (or 
in the direct process NN — > NNn) and absorption in the collision with a nucleon ttN — > R 
(or with two nucleons ttNN — > NN) is described by the l.h.s. of Eq . (|30|) . In BUU the 
test particle method is used to solve the transport equations. The pion spectral distribution 
function is represented by test particles in the following way: 

F n (r, t, k, Ml) = £ S(t - r< (t))(5(k - Ut))KMl ~ <,(*)) • (38) 



Following Refs. |], [16] we introduce a scalar off-shell potential acting on the i-th pion 
test particle: 

where p^(r,t) is the density of nucleons in the local r.f. of nuclear matter. t cr is the 
production time of the test particle. The potential (|39[) drives the off-shell pion back to its 
mass shell when it leaves the nucleus, i.e. we assume that 

M Wji {t)=m n + s^ i {r i {t),t) . (40) 

Therefore, Eqs. (|39|) J40|) provide the correct boundary condition for large times, i.e. when 
pions reach the vacuum: 

lim M„ At) = . (41) 

t — >oo ' 

The operator L ff^ s h e ii[F n ] in the l.h.s. of Eq.(^) is defined by its action on the distribution 
function in the test particle representation fl3~8|) as follows: 

d dH ff^ s heii,i d \ 



J of f— shell 



[F n ] = Yl \ 2M^i(t)s^ 



dMl dri dkj 

x 5(v - r,(t))*(k - k,(t))5(M 2 - Ml -{t)) , (42) 



where 



H ff shell, i = Jkf + {m n + s n)i ) 2 (43) 
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is the energy of the pion test particle i. From Eqs.(|30f) , (|42j) we get equations of motion of 
test particles in coordinate and momentum space: 



dH off 

—shell, i 



(44) 



Ok, 

(45) 



shell, i 



These two Hamiltonian equations together with Eq.(^) determine the evolution of the 



spectral distribution function (38) completely. 

Due to the explicit time dependence of the potential (^), which acts as an external 
field, the energy of the whole system is not conserved. However, in practice, the violation 
of the energy conservation due to the off-shell potential turns out to be very small (~ 0.2 
MeV/nucleon for the central Au+Au collision at 1 A GeV, i.e. ~ 0.1% of the total cm. 
kinetic energy per nucleon). This is caused by a relatively small number of pions with respect 
to the baryon number. 

Some comments are in order on the relation of the phenomenological approach of Refs. 
|9|, Kj applied in the present work with the more general test particle equations of motion 



derived in Refs. |10|, |21| on the basis of the Kadanoff-Baym equations. Indeed, Eqs.(p9|).([40D 
which determine the evolution of the pion mass and the Hamiltonian equations ([4*4]) , fl4li|) can 
be obtained from the formalism of Refs. [II], |2T|] by assuming that Im n c oc p^ and neglecting 
Ren c in the propagation of the pions. Earlier BUU studies [IT, [T2|, T3|] within the on-shell 



approximation have shown that the inclusion of the real part of the pion optical potential 
results in an enhanced soft pion yield. However, our attempted treatment of the off-shell 
pion propagation within the formalism of Refs. [I(], with the A-hole pion polarization 
function resulted in the appearance of the superluminous pions. We have, therefore, chosen 
the simplified Hamiltonian approach of Refs. || |16[ for the propagation of the off-shell pions 
which has no this problem, since the off-shell potential j is momentum-independent. 



IV. NUMERICAL RESULTS 



The pion spectral function has been implemented in the BUU model in the version of Ref. 
|i~6| . The temperature T and the baryon density p have been calculated locally in BUU and 
their values have been used to interpolate A n which was stored on the 4-dimensional grid in 
the space (M%, |k|, p,T). This allows one to save CPU time by avoiding on-line calculation 
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of A w . Masses of outgoing pions from resonance decays or from NN collisions have been 
sampled according to A w . The pions have then been propagated by using Eqs.(PD,(pE5D until 
they have been reabsorbed in 7rN — > R or ttNN — > NN processes or emitted to vacuum. The 
BUU calculations of the present work are performed using the soft momentum-dependent 
mean field (SM) with the incompressibility K = 220 MeV and employing the quenching 
of the resonance production in order to get the correct pion multiplicity Below, if 

the opposite is not stated explicitly, the vacuum A-width ([L4]) is used in the pion spectral 
function. 

Fig. |6| shows the time evolution of the pion mass distribution in the central collision 
Au+Au at 1.06 A GeV. At early times (t = 10 -j- 20 fm/c) a rather broad range in M n 
is populated by the resonance decays in dense nuclear matter, where the pion spectral 
function is broad. In order to separate the effect of the off-shell transport by Eqs.([44|).(|45|) 
we performed also the calculation putting L ff_ she u[F n ] = in Eq.(|30"D (dash-dotted lines). 
In this case pions emitted to vacuum do not lose their off-shellness and the M^-spectrum 
stays quite broad even at t = 40 fm/c when all violent processes are already finished and 
the system is already quite dilute (central baryon density ~ 0.02 fm~ 3 ). The full calculation 
(solid lines) results in a strong narrowing of the M^-spectrum at the late stage of the reaction: 
at t = 40 fm/c the most part of pions is already on the mass shell. Continuing the calculation 
until 80 fm/c (dotted line on lower right panel) produces a rather moderate variation of the 
M^-spectrum towards precise 6(M 1T — m K ). Thus, for the comparison with experimental data 
below we have stopped calculations at t = 40 fm/c. 

In Fig. we present the inclusive transverse momentum spectra at midrapidity for the 
same system Au+Au at 1.06 A GeV in comparison with the experimental data from Refs. 
|23| , |24| . The calculation with on-shell pions (dashed lines) describes reasonably well the high- 
p t part of the spectra, but underpredicts the data at low p t . Taking into account the pion 
off-shellness (dotted lines) increases the low-pt pion yields improving, thus, the agreement 
with the data. The low-pt pion enhancement is caused by the off-shell pion propagation. 
To check this, we switched off again the L ff- s h e u[F w ] term in Eq.(^) (dash-dotted lines). 
We see that, indeed without this term the \ow-p t enhancement disappears. The effect is 
caused by the low mass pions (M n < m n ) produced by the resonance decays at early stage 
of the collision (Fig. []). These pions experience the attractive off-shell potential (|3"5|), and, 
therefore, they are decelerated propagating out of the dense region of the nuclear system. 
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One may notice also in Fig. |7| a slightly enhanced yield of high-pt pions by the off-shell 
calculation with vacuum A-width. The agreement with data at large transverse momenta 
can be improved by using the A-spreading width (solid lines). This is caused by the reduced 
strength of the pion spectral function at large pion masses (c.f. Fig. [5]). 



V. SUMMARY AND OUTLOOK 



In this work we have generalized the BUU approach for the pion off-shellness in the nuclear 
medium. The A-hole model has been applied to calculate the pion spectral function which 
has been used, then, in BUU to choose randomly the outgoing pion mass in the resonance 
decay. The propagation of the produced off-shell pions has been described in a Hamiltonian 
approach - developed earlier in Refs. |)], |16[ - using the baryon density dependent off-shell 
scalar potential which acts on pions. This guarantees that a pion emitted to vacuum is 
always on the mass shell. 

As an observable effect in heavy-ion collisions at 1 A GeV, we have demonstrated the low 
p t enhancement of the pion yield due to the off-shell pion propagation, which is in agreement 
with earlier BUU calculations Ulli |i~3|l using the in-medium pion dispersion relation. We 
expect, that the pion off-shellness will also influence other observables, which are not covered 
in this work, e.g. the K + , K~ and p production via reactions ttN — > K + A, ttA — > K~N 
and 7T7T — > p respectively (study in progress). 

Another important aspect is the modification of the A-width due to coupling to an 
off-shell pion. We have studied this effect assuming the in-medium broadening of the A- 
resonance given by the phenomenological A-spreading potential ||18|| . Overall, a rather weak 
influence of the in-medium A-width on pion observables is found. We have to remark, 
however, that we included the in-medium A-width in the pion spectral function only, not in 
BUU. The explicit in-medium modification of the A-width in BUU has been studied earlier 
in Refs. [|11], |25|, where a quite weak effect on pion observables in heavy- ion collisions and 
on 7-absorption cross section on nuclei was found. Finally, it is straightforward to generalize 
the A-hole model including higher baryonic resonances which will then also contribute to 
the pion spectral function. 



14 



Acknowledgments 

Stimulating discussions with W. Cassing, S. Leupold, M. Post and V. Shklyar are grate- 
fully acknowledged. The authors are especially grateful to Stefan Leupold for the careful 
reading of the manuscript, useful advices and comments. 



15 



[1] W. Ehehalt, W. Cassing, A. Engel, U. Mosel, and Gy. Wolf, Phys. Rev. C 47, R2467 (1993). 

[2] S.A. Bass, C. Hartnack, H. Stocker, and W. Greiner, Phys. Rev. C 51, 3343 (1995). 

[3] A.B. Migdal, O.A. Markin, I.N. Mishustin, Sov. Phys. JETP 39, 212 (1974). 

[4] G.E. Brown and W. Weise, Phys. Rep. 22, 279 (1975). 

[5] T. Ericson and W. Weise, "Pions and Nuclei", Clarendon Press, Oxford, 1988. 

[6] PA. Henning, H. Umezawa, Nucl. Phys. A 571, 617 (1994). 

[7] R. Rapp, J. Wambach, Nucl. Phys. A 573, 626 (1994). 

[8] G.F. Bertsch, P. Danielewicz, Phys. Lett. B 367, 55 (1996). 

[9] M. Effenberger and U. Mosel, Phys. Rev. C 60, 051901 (1999). 

[10] W. Cassing, S. Juchem, Nucl. Phys. A 665, 377 (2000); Nucl. Phys. A 672, 417 (2000). 

[11] J. Helgesson and J. Randrup, Phys. Lett. B 439, 243 (1998). 

[12] L. Xiong, CM. Ko, and V. Koch, Phys. Rev. C 47, 788 (1993). 

[13] C. Fuchs, L. Sehn, E. Lehmann, J. Zipprich, and Amand Faessler, Phys. Rev. C 55, 411 
(1997). 

[14] M.J. Dekker, P.J. Brussaard, J.A. Tjon, Phys. Rev. C 49, 2650 (1994). 

[15] M. Post, S. Leupold and U. Mosel, Nucl. Phys. A 689, 753 (2001). 

[16] M. Effenberger, E.L. Bratkovskaya, and U. Mosel, Phys. Rev. C 60, 44614 (1999). 

[17] S. Teis, W. Cassing, M. Effenberger, A. Hombach, U. Mosel and Gy. Wolf, Z. Phys. A 356, 
421 (1997). 

[18] M. Hirata, J.H. Koch, F. Lenz and E.J. Moniz, Ann. of Phys. 120, 205 (1979). 

[19] E. Oset and L.L. Salcedo, Nucl. Phys. A 468, 631 (1987). 

[20] CM. Ko, L.H. Xia and P.J. Siemens, Phys. Lett. B 231, 16 (1989). 

[21] S. Leupold, Nucl. Phys. A 672, 475 (2000). 

[22] A.B. Larionov, W. Cassing, S. Leupold, U. Mosel, Nucl. Phys. A 696, 747 (2001). 

[23] D. Pelte et al., Z. Phys. A 357, 215 (1997). 

[24] O. Schwalb et al., Phys. Lett. B 321, 20 (1994). 

[25] M. Effenberger, A. Hombach, S. Teis, U. Mosel, Nucl. Phys. A 613, 353 (1997). 



16 




FIG. 1: Feynman graph representing — iH(k) with n(fc) being the pion polarization function, k 
and p are the pion and nucleon four-momenta respectively. 
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FIG. 2: Real (upper panel) and imaginary (lower panel) parts of the pion optical potential as 
functions of the pion momentum at various densities and temperatures. 



A n [GeV" 2 ] 




FIG. 3: Pion spectral function in nuclear matter at the density p = po and temperature T = 5 
MeV shown vs mass squared and momentum. 




FIG. 4: Peak positions of the pion spectral function A n in the (|k|, w)-plane for nuclear matter 
at the density p = po and temperature T = 5 MeV - full circles. Free pion dispersion relation 
E w (k) is shown by the solid line. Dashed line shows the A-hole energy 2£jvA(k). Solution et of the 
dispersion relation (f2^) is depicted by dash-dotted line. 
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FIG. 5: Pion spectral function vs pion invariant mass squared at fixed pion momentum |k| = 0.25 
GeV/c (upper panel) and |k| = 0.55 GeV/c (lower panel) in nuclear matter at p = 2po and 
T = 70 MeV. Dashed and solid lines present the calculations with vacuum and in-medium A-width 
respectively. 
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FIG. 6: Time evolution of the pion mass distribution in the central Au + Au collision at 1.06 A 
GeV. Solid and dash-dotted lines show the results of calculations with and without the off-shell 
term L ff- S h e ii [F n ] in pion propagation respectively. Dotted line in the lower right panel shows 
the mass distribution at t = 80 fm/c for the calculation with off-shell term. 
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FIG. 7: Pion transverse momentum spectra at midrapidity for Au+Au collisions at 1.06 A GeV 



in comparison with experimental data from Refs. 24]. Dashed lines show the calculation with 
on-shell pions. Dotted and dash-dotted lines present the off-shell calculations using vacuum A- 
width with and without L ff_ s f ie n[F 7T ] term respectively. Solid lines show the off-shell calculation 
with in-medium A-width. 



